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^ A new limit is presented on the axion-like monopole-dipole P,T- 

non-invariant interaction in a range (10 — 1) cm. The spin- 

dependent nucleon-nucleon potential between neutrons and nu- 
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^ 1 Introduction 

Hypothetical pseudoscalar particle - axion offers a window to probe very small 
coupling and very high energy scales [1]. 

Axion, according to later modifications of the primary model [1], has mass in 
a very large range: (10^^ < rua < 10^) eV. Current algebra technique is used 
to relate the masses and coupling constants of the axion and neutral pion: = 
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{f^rriT^/ f a)^/z/{z + 1), where z = mu/nid = 0.56, f^^ « 93MeV, = 135 MeV, 
so that nio, ~ (0.6 x lO^'^ GeV/fa) meV. Here is the scale of Peccei-Quinn 
symmetry breaking. The axion coupling to fermions can be in general represented 
as Qaff = Cfirif/ fa, where Cf is the model dependent factor [2,3]. 

Early reactor, beam-dump, weak decays, and nuclear transition experiments 
have placed lower limits on the axion mass. Stringent limits, especially from the 
lower side of the axion mass range, have been set on the existence of axion using 
astrophysical and cosmological arguments [2,4]. These more recent constraints 
limit the axion mass to (10~^ < nia < 10^) eV with correspondingly very small 
coupling constants to quarks and photon [2-A]. Although these limits are more 
stringent than can be reached in laboratory experiments, it is of interest to try to 
constrain the axion as much as possible using laboratory means. Interpretation of 
laboratory experiments depend on less number of assumptions than the constraints 
inferred from astrophysical and cosmological observations and calculations. The 
laboratory experiments performed or proposed so far are rather diverse and employ 
a variety of detection techniques. Axion is still one of the candidates for the cold 
dark matter of the Universe [5]. Some recent reviews are [2, 3]. 



2 IVIonopole-dipole interaction potential 

Axions mediate a P- and T-odd monopole-dipole interaction potential between 
spin and matter (polarized and unpolarized nucleons) [6]: 

where Qs and gp are dimensionless coupling constants of the scalar and pseu- 
doscalar vertices (unpolarized and polarized particles), m„ is the nucleon mass 
at the polarized vertex, cr is vector of the Pauli matrices related to the nucleon 
spin, r is the distance between the nucleons, A = h/iriaC is the range of the force, 
TUa is the axion mass, and n = r/r is the unit vector. 
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The microscopic potential ([T]) between two nucleons creates a macroscopic po- 
tential between nucleon and matter. If the matter is represented by a layer of 
thickness d, the neutron interaction with it is described by the potential 

V{x) = {a . n)g,gj^e~^'\l - e"^/^) = V,{ct ■ n)e~^'\l - e"'^/^), (2) 

where x is the distance to the interface along the normal (unit vector n), is the 
nucleon density in the layer, and in the last equality we defined 

We will deal here only with ultra-cold neutrons (UCN) or total reflection of higher 
energy neutrons, and because of that we do not consider how the potential changes, 
when a neutron enters the matter. 

Several laboratory searches provided constraints on axion-like coupling in the 
macroscopic range A > 0.1 cm [2,7]. The limit in the A-range (10^ — 1) cm was 
established in the Stern-Gerlach type experiment in which UCN were transmitted 
through a slit between a horizontal mirror and absorber [8]. The obtained limit 
for the value QsQp was ~ 10^^ at A = 10^ cm, that corresponds to the value 
of the monopole-dipole interaction potential at the surface of the mirror Vq ~ 
10^ neV. This value is equivalent to an effective magnetic interaction —fiB of 
the neutron with magnetic field B ~ 0.2 G. It is estimated that in future ultra-cold 
neutron Stern-Gerlach type experiments sensitivity will be improved by orders of 
magnitude [9]. A better sensitivity is also expected in a proposed experiment on 
the ultra-cold neutron magnetic resonance frequency shift [10]. 

We consider here what limits on QsQp in the range A (10^ — 1) cm can be 
extracted from depolarization of UCN in storage traps. Depolarization can be 
expected because the particle spin interaction with axion field Uax = (cr ■ n)V is 
similar to magnetic interaction Umagn = ' B) and corresponding pseudo- 

magnetic field in general is not collinear to the neutron polarization. 

Depolarization was already considered earlier in the paper [11], however it 
was estimated there semiclassically as the neutron spin rotation in the interaction 



region A in vicinity of the reflecting waU, which is not sufficiently rigorous. Here 
we calculate the depolarization probability at a single collision with the wall with 
distorted wave Bom approximation (DWBA) perturbation theory. 

3 Constraints from the ultra-cold neutron depolarization 

In fact, depolarization of UCN in traps can be attributed to many factors. In be- 
tween them are inhomogeneities of the internal magnetic field and presence of 
magnetic impurities on the walls. We use here the most conservative estimate at- 
tributing all the neutron depolarization to the hypothetical axion-like interaction. 

Let's consider a nonmagnetic semi-infinitely thick wall with optical potential u 
at X > 0, external homogeneous field B parallel to z-axis, and the axion pseudo 
magnetic field h{x) = 6oexp(x/A) parallel to x-axis. We want to calculate spin- 
flip reflectivity of a neutron initial polarization along the B field, taking the ax- 
ion field as a perturbation. For that we need to solve one dimensional stationary 
Schrodinger equation, which can be represented in the form 

[(f/dx'^ + - Ba, - h{x)a^e{x < 0) - ue{x > 0)]|^(a;)) = 0, (4) 

where a^^z are the well known Pauli matrices, Q{x) is a step function equal to 
unity, when inequality in its argument is satisfied, and to zero otherwise, k is 
normal component in vacuum of the wave vector of the incident particle, and 
the magnetic fields include factor 2m\ii\/h^, which contains neutron mass m and 
magnetic moment ii. To find a spinor solution |^(a:)) of (|4]) we need to define an 
incident wave \iIjo{x)). In general we can define it as 

lijjoix)) = exp{ikx)\^) , (5) 

^ . 

where k = y/c^ — Baz, and |^) is an arbitrary spin state, which is a superposi- 
tion |,^) = a\^u) + P\^d) of states 1^^,^) that are eigen spinors of the matrix a^: 
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A non-perturbed solution of (|4]) is 

l^o(x)) = |e(x < 0)[e*^" + e-*^^r(A;)] + e{x > 0)i{k)e''^'A |0, (6) 



where k = -s/k^ — Bg^ — u, reflection, r, and transmission, t, matrices deter- 
mined (see, for instance [12]) from matching conditions at the interface are 

\ ki /c ^, , 

m = . — ^, t{k) = ^ — ^. (7) 

k^k k+k 

It is seen that the spinor ^ can be represented as |^o(a^)) = ^o(a^)IO' where 
^q{x) is 2 X 2 matrix 

#0(2^) = < 0)[exp(zfcx) + exp(-zfcx)r(/c)] + 6(x > 0) exp(zfc'a:)f(A;). (8) 

This matrix satisfies the non-perturbed Schrodinger equation 

[d^/dx^ + /c^ - Ba, - ue{x > Qi)]^q{x) = 0, (9) 

and is diagonal one, which means that its non-diagonal matrix elements 

^Q{x)du = {U^Q{x)\iu). ^Q{x)ud = {iu\^Q{x)\id) (10) 

are zero. 

The perturbation h{x)(7x changes and the change \5^{x)) = 5'^{x)\C), 
according to perturbation theory is representable as 

5^{x) = I G{x,x')b{x')a^'^o{x')dx', (11) 

where the matrix Green function, G, in the DWBA approach is a causal solution 
of the inhomogeneous Schrodinger equation: 

[d^/dx^ + k^ - Ba, - wG(x > 0)]G{x, x) = 16{x - x). (12) 

Here / in the right hand side is the unit matrix. Solution of this equation is con- 
structed with the help of two linearly independent solutions ^1^2(2^) of (M): 



G{x, x) = w-\e{x > x)^i{x)^2{x) + e{x < X )*2(a:)*i(x )], (13) 



where w is their Wronskian 



(14) 



and prime means derivative over x: ^ {x) = d^{x)/dx. 

For ^i(a:) we can take solution (8): ^i(x) = 4fo{x), and for linear indepen- 
dent solution ^2{x) we can take 

^2(2:) = 0{x < 0) exp{—ikx)t {k) + 6(x > 0)[exp(— ifc x) + exp{ikx)f'{k)], 

(15) 

where matching conditions satisfy for 



k -k 



r\k) = -. ^ = -r(/c), t{k) = -. 



2k 



k + k 

The Wronskian of ^1^2(2^) is equal to 

w = 2ikt{k). 



^1 • 



k + k 



(16) 



(17) 



Substitution of all these matrices into (11) gives 6'^{x) = exp{—ikx)R, where 



R 



-00 



2ik 



(18) 



The searched amplitude of spin flip reflection is 



Rdu = {^d\R\i 



U 5 



(19) 



Substitution of ( |T8| ) gives 

^ dx' 



R 



du 



[e'^''^' + e-'^"^'rd{k)]h{x)[e"'-^' + e-*""^V„(A;)], (20) 



J\ \Jkux' 



2'lkn 



where ku4 = Vk^T B, 



^r' , k,,, = VW^, /c;,, = VP - u T i^. (21) 



Substitution of 6(a:) = 60 exip{qx), where q = 1/A, and integration over x' in (20) 
gives 



Rdu — 



1 



+ 



rurd 



2ikd \q + i{ku + kd) q- i{ku + kd 



^ rd ^ ru 

q + i{ku -kd) q- i{ku - kd] 

bo (q{l + rurd) -i{ku + kd){l-rurd] 



2ikd V 9^ + {K + kd^ 

q{rd + ru) - ijku - kd){rd - r^) \ 

In the case of total reflection and not large external field {Umagn = ' B) <C 
En) we can approximate ku + kd ~ 2A;, ku — kd ~ —B/k and Vu ^ Vd = 



exp{—2i(j)), where = arccos(/c/y^). In this approximation the Eq. (22) is 
reduced to 

o _ ,-2z0^ ( gcos(20) + 2/csin(20) g \ 

As we are interested in the interaction range satisfying to /cA 1 (typical 
UCN k ~ 10^ cm~\ A > 10 cm) and not too strong external magnetic fields 
(<~ 500 G), the first term in Eq. (23) can be neglected comparing to the second 
one: 

and the spin-flip reflectivity from the wall, represented in dimensional units, be- 
comes 

- , iVj < fx (1 - e-^'^f 
w — \±tdu\ — r2 5 

where cjq = ")nB is the neutron spin Larmor frequency in the external field B, 
7„ = 1.83 X 10"^ s VG - the gyromagnetic ratio for the neutron, < i;^ > is the av- 
eraged over the UCN spectrum normal to the surface neutron velocity component. 

At cjo ^< ^± > /A we have the expression coinciding with the one obtained 
for the quasiclassical case and derived in the Ref. [13] with loq^ behavior of the 
depolarization probability. Typically < vj^ >~ 300 cm/s, and for A = 10^ cm 
the quasiclassical case is valid at 5 > 100 G. 

At a weak guiding magnetic field ~ 10 ^ G the quasiclassical approach is valid 
only for the interaction range Kx vj^> /cjo^l cm. 



Substituting (2) into (25) we obtain the expression for gggp. 

where (3 is the experimentally measured UCN depolarization probability per one 
reflection from the wall of the storage cavity. 

There are two published experimental data on the ultra-cold neutron depolar- 
ization: [14] and [15], in which special experiments are described to measure this 
value. In both publications, for a variety of materials, the measured values of 
the neutron depolarization probability per one neutron collision with the walls of 
storage cavity were around P ~ 10^. The (3 ~ 10^ was measured in [15] for 
the diamond like carbon foils (DLC). The magnetic fields in the storage chambers 
were partly due to stray fields from strong magnets used for the polarization of the 
incident neutrons [14, 15] and partly were formed by special magnets [15]. In [14] 
magnetic field was estimated [16] to be near 50 G , and in the experiment [15] it 
was reported to be ~55 G. In both cases, this rather large external magnetic field 
suppressed effect of the additional hypothetical spin-dependent interaction on de- 
polarization of the ultra-cold neutrons in traps, and therefore decreased sensitivity 
of these measurements to establishing constraints on the axion-like interaction. 

A better constraints can be obtained from the measurements of the ultra-cold 
neutron depolarization in traps at lower magnetic field in the experiments on 
the search for the neutron electric dipole moment (EDM) [17] and [18]. There 
the ultra-cold neutrons preliminary polarized by transmission through magnetized 
ferromagnetic foil were stored in a cylindrical bottle permeated by magnetic and 
electric fields. The magnetic field was applied parallel to the axis of the bottle, 
and its value in these experiments was very low: Bz = 0.02 G in [17], and Bz = 
0.01 G in [18]. The change of the magnetic resonance frequency was sought for 
at the reverse of the electric field direction. After filling the bottle with ultra- 
cold neutrons and closing the neutron valve, the 7r/2 Ramsey pulse was applied, 
which turned neutron spins perpendicular to the magnetic field. The neutrons were 
allowed to precess about magnetic field for 130 s [18], after which the second 7r/2 
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Ramsey pulse was applied, and neutrons in the appropriate spin state passed back 
through the polarizing foil to the neutron detector. 

Depolarization of neutrons at reflections from the walls of the storage cav- 
ity in presence of a gradient of a spin-dependent potential decreases contrast 
of the neutron magnetic resonance curves. Probability of the neutron depolar- 
ization was not measured directly in these experiments, but from the reported 
very good magnetic resonance curves it can be concluded that the UCN depo- 
larization probability at a single collision with the walls was not higher than 
in [14, 15]. According to [19] the neutron depolarization time in the EDM ex- 
periment [18] can be estimated to be not less than Tdep ~ 800 s. From this figure 
we can estimate the depolarization probability per single collision with the walls: 
P = l/{rdepv) = 18/(800 X 500) ~ 4 x 10'^, where / is the mean free path 
between two consecutive collisions, which can be found from dimensions of the 
neutron bottle: in the experiment [18] it was a cylinder of the internal diameter 
about 44 cm and height 15 cm. We constrain the monopole-dipole interaction in 
two assumptions for the depolarization probability in [18]: one - /3 = 4 x 10~^, 
and another one - corresponding to the best value obtained in the experiments of 
the PSI group [15] - /? = 10"^ 

With above f3 we can draw the limiting curves for the parameters of the monopole- 
dipole coupling of the axion field, which is shown in Fig. 1 together with results 
from other publications. 

The ultra-cold neutron depolarization data may be used also to set limits on 
the monopole-dipole coupling between neutrons and electrons of the walls of the 
storage chambers. However, because density of the electrons in the medium is 
approximately two times lower than the density of nucleons, the constraints are 
respectively two times less strong. 
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X, cm 

Figure 1: Constraints on the axion monopole-dipole coupling strength gsQp and effective range: 1 
and 2 - constraints for the value of coupling constant of nucleon and electron gsQp from Refs. [20] 
and [7], respectively; 3 - from the neutron gravity experiment of Ref. [8], 4 - from Ref. [11], 
5 - from spin relaxation of ^He, Ref. [21], 6 - this work at the UCN depolarization probability 
f3 = 10"^ and magnetic field = 50G [15], 7 - the same, but B, = 0.01 G [18]; 8 - the same, 
but Bz = 0.01 G, P = 4: X 10"^ [18, 19]. It was assumed in both cases of the ultra-cold neutron 
storage, that d = 1 cm, N ^ 2 x 10^^ cm"^, < >= 300 cm/s. 
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